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Abstract 

We study the combined effects of electrostatic and 
hydrodynamic interactions (HI) on the short-time 
dynamics of charge-stabilized colloidal spheres. For 
this purpose, we calculate the translational and the 
rotational self-diffusion coefficients, £>* and D^, as 
function of volume fraction </> for various values of 
the effective particle charge Z and various concen- 
trations Ug of added 1-1 electrolyte. 

Our results show that the self-diffusion coeffi- 
cients in dcionized suspensions are less affected by 
HI than in suspensions with added electrolyte. For 
very large Ug, we recover the well-known results for 
hard spheres, i.e. a linear (/)-dependence of D\ and 
Z)g at small (j). In contrast, for deionized charged 
suspensions at small 0, we observe the interesting 
non- linear scaling properties D\ cx 1 — at(j)^^^ and 
Dg oc 1 — a^f/)^. The coefficients at and are found 
to be nearly independent of Z. The qualitative dif- 
ferences between the dynamics of charged and un- 
charged particles can be well explained in terms of 
an effective hard sphere (EHS) model. 

1 Introduction 

Since several years, the effect of HI on the short-time 
self-diffusion coefficients of hard sphere suspensions 
has been investigated in detail by various authors 
[|[]5|. For the calculation of the first and second 
virial coefficients of Z)* and in an expansion in 
terms of the volume fraction 0, both the influence of 
two-body and three-body HI was taken into account. 
At small <j), the currently established results for the 



normalized diffusion coefficients ff* and are given 
by 0,0 
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and by @| 



Pl 



= l- 0.630(/)- 0.67(/)2-f 0((/|3), (2) 



respectively. Here, Dq and Dq are the Stokesian dif- 
fusion coefficients for a colloidal sphere of radius a 
dispersed in a solvent of viscosity rj. 

The possibility to express if* and iJJ in terms of a 
power series in (p arises from the fact that hard sphere 
suspensions at small (p can be considered as dilute 
both with respect to the particle hydrodynamics and 
to the microstructure. For charge-stabilized suspen- 
sions, however, this is not possible in general 1^,0]. 
Especially deionized, i.e. salt-free suspensions ex- 
hibite pronounced spatial correlations even at very 
small (f), so that these systems are diluted only as 
far as the HI is concerned. The corresponding ra- 
dial distribution function g{r) has a pronounced ip- 
dependence, a well developed first maximum, and it 
shows a correlation hole, i.e. a sperical region with 
zero probability for finding another particle, which 
usually extends over several particle diameters [^|^ . 
In contrast, the g{r) of hard spheres is nearly a 
unit step function g{r) 0(r — 2a) for < 0.05. 
Therefore, the calculation of iJ* and -ffj at small 
(j) is more demanding for charged suspensions than 
for hard spheres, because for the charged particles 
it is necessary to use distributions functions gener- 
ated from computer simulations or integral equation 
methods 

We will show subsequently, that it is essentially 
the presence of the correlation hole for charged sus- 
pensions, which causes large and interesting differ- 
ences in the short-time diffusion of charged and un- 
charged suspensions. 



1 



2 Calculation of and if^ 

s s 

In the following, we shortly summarize the main ex- 
pressions needed to calculate _ff * and _ff J for charge- 
stabilized suspensions. A more detailed description 
of the method used by us for the calculation of short- 
time diffusion coefficients is given in Refs. @,|]. 

As shown in Refs. j^,^, both and iJJ can be 
measured using depolarized dynamic light scatter- 
ing (DDLS) from suspensions of optically anisotropic 
colloidal spheres. On the time scales, which are ac- 
cessable by DDLS, the theoretical expression for i?* 
is given by H 

Hl^^^{TTD{\ir^)). (3) 

The corresponding expression for iJJ is obtained 
from eq. (^) by simply replacing the superscript 
't' by v. The hydrodynamic diffusivity tensors 
D*i and relate the force/torque exerted by the 
solvent on an arbitrary particle 1 with its transla- 
tional/angular velocity |3[0. Due to the many- 
body character of HI, both tensors depend on the 
particle configuration — (ri,...rjv) of all N 
interacting particles, and in principle the full A^- 
particle distribution function is needed to perform 
the ensemble average (...}. TrJ)"]^ denotes the sum 
over the diagonal elements of . 

For an appropriate evaluation of eq. (^, we use 
a rooted cluster expansion || , which leads to a 
"hydrodynamic virial expansion" of iJ*: 

Hl^l + Hl,<P + Hl,q^' + Oiq^'). (4) 

Here, the coefficient iJ*]^ is given by an integral over 
the product of g{r) with a translational hydrody- 
namic mobility function, which depends only on the 
distance r of two spheres [|l|,|^,|7[|| . The second coeffi- 
cient accounts for three-body HI. For evaluating 
iJ*2 1 o'^G needs therefore an expression for the static 
triplett correlation function ^^'^•'(r, r') which appears 
as part of the integrand of a three-fold integral. 

A similar analysis is used to calculate lead- 
ing to results which involve now rotational hydro- 
dynamic two-body and three-body mobility func- 
tions 

The results for Hi and depicted in eq. (|l|) 
and (||) were derived from eq. @ by using in iJ*^ 
and the g{r) of hard spheres evaluated up to lin- 
ear order in </>, whereas the vanishing density form 
of g^^\r,r') was used in calculating the coefficients 
and -ffj2- I'^ these results, exact two-body HI is 
accounted for Hl^ and iJJ^, whereas only the leading 
long-distance contribution to the three-body mobil- 
ity functions was used for the calculation of 7J*2 

For charged suspensions, however, it is not pos- 
sible to use in eq. (g) low-order virial expres- 
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Figure 1 : Normalized short-time translational diffusion 
coefficient Hi for a deionized charged suspension with 
Z = 200, a = 45nm, T = 294K, and e = 87.0. Also 
shown is the result for hard spheres according to eq. (|l|). 

sions of the two-body and three-body static distri- 
bution functions. In this study, we use instead re- 
sults for g{r), which are obtained from the rescaled 
mean spherical approximation (RMSA), as applied 
to the one-component macrofluid model of charge- 
stabilized colloidal suspensions ||]. The effective 
pair potential u{r) acting between two particles is 
modelled by the repulsive part of the famous DLVO- 
potential, i.e. (3u(r) — ATexp [— K(r — 2a)] -, for r > 
2a. Here, /3 = {kBT)-\ K = Z^{LB/a){l + na)'^, 
Lb = /3e^/e, and e denotes the dielectric constant 
of the solvent. The screening parameter k is given 
by = Lb [3|Z|0/a^ -I- Snus] , where Hs is the con- 
centration of added 1-1 electrolyte, and the counte- 
rions are assumed to be monovalent Moreover, 
we use Kirkwood's superposition approximation for 
g^^\r,r'), inserting again the RMSA-(7(r). Further 
details concerning the numerical calculation of iJ^i, 
ffJi, HI2, and are given in Refs. |,|. 

3 Results and discussion 

We focus first on the short-time diffusion coefficients 
of deionized charged suspensions, i.e. where Us = 0. 
Our results for Hi and i?J are shown in figs. |^ and 
^ The used system parameters are typical for sys - 
tems which have been under experimental study |ll|| . 

Obviously, the effect of HI on the self-diffusion co- 
efficients is less pronounced for charged suspensions 
than for hard spheres at the same (f>. Furthermore, 
we find a quite different volume fraction dependence 
of Hi and i/J for charged and uncharged particles. 
Whereas for hard spheres the 0-depcndencc of Hi 
and iJJ is linear at small (cf. eqs. (|^) and (||)), 
we obtain from a least-square fit of our numerical 
results (shown as crosses in figs. and 1) the fol- 
lowing results for deionized charged suspensions for 
< </> < 0.05 g 

Hi = l-at0i-3O, at = 2.59, (5) 
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Figure 2: Normalized short-time rotational diffusion 
coefficient ffj for a deionized suspension with system 
parameters as in fig. |l|. Also displayed is the result for 
hard spheres according to eq. (Q). 
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/1. 99 



1.28. 



(6) 



The coefficients at and are found to be nearly 
independent of the effective particle charge when 
Z > 200 @,|. Note from fig. |, that eq. (§ con- 
stitutes the best fit function for ffJ((/>) even in the 
extended interval O<0<O.15. In case of Hi how- 
ever, the parametric form if* = at(j)^ provides no 
good fit for values of extending beyond 0.05 ||]. 

There is a simple physical explanation for the 
weaker influence of HI on the self-diffusion coeffi- 
cients of charged suspensions as compared to un- 
charged ones. As already mentioned, the g{r) of 
deionized suspensions displayes a pronounced cor- 
relation hole, resulting form the strong electrostatic 
interparticle repulsion. Consequently, the hydrody- 
namic coupling between the translational or rota- 
tional motions of two spheres becomes rather small, 
thus giving rise to the observed weak influence of HI. 
Unlike charged particles, the influence particularly of 
the short-range part of HI is rather strong for hard 
sphere suspensions at small 4>- This is due to the 
large probability of finding hard sphere particles at 
contact or close to the contact distance r = 2a. 

Along this type of arguments, it is also possible to 
explain the differences of Hl{(j)) and Hl{4)) in deion- 
ized suspensions and suspensions with nonvanishing 
Hs- We only show here the results of our calcula- 
tions of H'g{(f>) for example. From these results in 
fig. H, we notice that iJJ becomes more and more 
affected by HI when rig is increased. For very large 
Us, HI of charged particles approches the result for 
hard sphere suspensions, obtained semianalytically 
in Ref. § (cf. eq. (|)). 

This finding is easily explained by noticing, that 
the extension of the correlation hole decreases with 
increasing n^, leading to a stronger hydrodynamic 
coupling of the particles. Upon addition of elec- 
trolyte, the electrostatic repulsion of the particles 
becomes more and more screened and short-ranged. 



Figure 3: Volume fraction dependence of HI for various 
amounts of added 1-1 electrolyte, as indicated in the 
figure. All other system parameters as in fig.lll. 



resulting in a pure hard-core repulsion for oo 
[^|-^. Therefore, the microstructure of the suspen- 
sion gradually transforms to that of hard spheres, 
with HI approching the parametric form given in 
eq. (|2|). We further note, that the radial distri- 
bution function g{r) corresponding to fig. ^ exhib- 
ites a small correlation hole even for — limM. 
This leads to the small differences of our results for 
ris — limM and for hard spheres in fig. ^. 

We mention, that our results for Hi show simi- 
lar trends, i.e. a gradual transformation of the 0- 
dependence of Hi from eq. (||) to eq. (0) with in- 
creasing ns ||]. 

In the remainder of this article, we focus on the 
qualitatively different ^-dependencies of Hi and HI 
found in case of deionized charged and uncharged 
suspensions. For an intuitive physical explanation, 
we use an effective hard sphere model (EHS model) 
[^0, describing the actual g{r) as a unit step func- 
tion gEHsi"!^) = ©(j' ^ '^CLEHs)- The EHS radius 
o-EHS > 0, accounts in a crude fashion for the corre- 
lation hole, observed in the actual g{r). We identify 
'2aEHS = fm, where Vm is the position of the first 
maximum of g{r). It is now crucial to notice that 
shows an interesting scaling property when rig = 0. 
Due to the strong electrostatic repulsion, has the 
same ^-dependence as the average geometrical dis- 
tance f between two spheres. Hence 
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f = a^47r/30"3. 



(7) 



Using the approximation gEHsi^) of 5('')i it is easy 
to calculate the coefficients Hl-^ and H^-^ in an ap- 
proximative way. By using far-field expansions of the 
hydrodynamic two-body mobility functions [f7|,p|,pX[|, 
one obtains the following results from the leading 
terms of these expansions 



Hi, = 
Hli = 



15 / a 
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(8) 
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This leads together with eq. (|^) and to the ex- 
pressions 

Hi = 1 - + 0{<j>^), A* > 0, (10) 

H: = 1- ^^^^^(^l)^ A'->0, (11) 

with exponents which are in good agreement with 
our numerical findings given in figs. |and| (cf. eqs. 
(i and d)). 

Therefore we have shown by a simple analytic 
calculation based on the EHS model, that the ob- 
served differences in the functional forms of Hl{(t)) 
and Hl{(j)) between charged and uncharged suspen- 
sions are mainly caused by the leading terms of the 
hydrodynamic two-body mobility functions in com- 
bination with the scaling property r,„ oc (j)^^^^ ^ valid 
for deionized suspensions. The higher order terms in 
the hydrodynamic far-field expansions only give rise 
to minor corrections to the observed scaling proper- 
ties depicted in eq. (||) and (||). These terms be- 
come increasingly important for larger volume frac- 
tions > 0.05 (cf. eqs. ^ and (|ll|) in the EHS 
model). 

When electrolyte is added to the suspension, eq. 
becomes invalid because of the enhanced screen- 
ing of the direct particle interactions. This causes 
a change in the functional behaviour of Hl{(t)) and 
{(/)), as can be seen both from the EHS model 
and from our numerical results (cf. fig. ^ in case of 
HI) 01. 

Using the EHS model, it is also possible the moti- 
vate the nearly Z-independence of at and in eqs. 
(||) and (||). Since r,„ is nearly independent of Z for 
Z > 200, the EHS model predicts charge indepen- 
dent results for the short-time diffusion coefficients 
of deionized suspensions, in agreement with our nu- 
merical results. 

We mention, that it is also possible to deal with 
and within the EHS model, giving further 
insight in the volume fraction dependence of the dif- 
fusion coefficients of deionized suspensions 0,1). It is 
then possible to explain qualitatively the surprising 
fact that -ffj(0) is well parametrized up to (/i = 0.15 
by the functional form ifj = 1 — a^^^, obtained in 
the EHS model by using only the leading term in the 
far-field expansion of the rotational two-body mobil- 
ity functions 0,1). 

4 Conclusion 

We have presented calculations of the translational 
and rotational short-time self-diffusion coefficients 
for charge-stabilized suspensions. The self-diffusion 
coefficients of charged suspensions are less affected 
by hydrodynamic interactions than the correspond- 
ing coefficients of hard spheres. As a major re- 
sult we have found substantially different volume 



fraction dependencies of i?* and iJJ for (deionized) 
charged and uncharged suspensions. The observed 
differences are well explained in terms of an effective 
hard sphere model by observing the big differences 
in the microstructure of suspensions of charged and 
uncharged particles. 

We note finally that recent DDLS measurements of 
-ffj in deionized suspensions of charged fiuorinated 
polymer particles compare favourably with our re- 
sults in eq. (^) jll). On the other hand, to our 
knowledge, no experimental data of iJ* for deion- 
ized charge-stabilized suspensions are accessable so 
far. We further point out that the interesting qual- 
itative differences between charge-stabilized suspen- 
sions and hard spheres exist also with respect to sed- 
imentation [HI and long-time self-diffusion [p^ . 
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